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Abstract: A combined time-domain fluorescence and hemoglobin diffuse 
optical tomography (DOT) system and the image reconstruction methods 
are proposed for enhancing the reliability of breast-dedicated optical 
measurement. The system equipped with two pulsed laser diodes at 
wavelengths of 780 nm and 830 nm that are specific to the peak excitation 
and emission of the FDA-approved ICG agent, and works with a 4-channel 
time-correlated single photon counting device to acquire the time-resolved 
distributions of the light re-emissions at 32 boundary sites of tissues in a 
tandem serial-to -parallel mode. The simultaneous reconstruction of the two 
optical (absorption and scattering) and two fluorescent (yield and lifetime) 
properties are achieved with the respective featured-data algorithms based 
on the generalized pulse spectrum technique. The performances of the 
methodology are experimentally assessed on breast-mimicking phantoms 
for hemoglobin- and fluorescence-DOT alone, as well as for fluorescence- 
guided hemoglobin-DOT. The results demonstrate the efficacy of improving 
the accuracy of hemoglobin-DOT based on a priori fluorescence 
localization. 
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1. Introduction 

Breast tumor diagnosis is widely regarded as one of the most applicable areas of diffuse 
optical tomography (DOT), which aims at reconstructing blood concentration and 
oxygenation images of the breast by spatially and spectrally resolving its tissue optical 
properties [1—5]. This endogenous-contrast-based imaging modality, referred to as 
hemoglobin-DOT, has potentials to compensate the limitations of the established modalities, 
such as X-ray mammography, ultrasound (US) and magnetic resonance imaging (MRI), in the 
safety, sensitivity and specificity, particularly for screening young women under 35 years of 
age who have dense breasts [1,4,5], and allows for improved discrimination between 
malignant and benign lesions [3,5-7]. Breast hemoglobin-DOT has been widely explored with 
three distinct technologies, following the same evolutional track as the universal diffuse light 
imaging: continuous-wave [8], frequency-domain [9], and time-domain (TD) [6,10]. Many 
groups have reported high total hemoglobin concentration in malignant tumors as compared to 
normal tissues - an expectable contrasting mechanism due to the sustained angiogenesis 
occurring in tumors [3,11,12]. Despite of unapparent manifestation in hemoglobin-DOT, 
some groups have observed a decrease of the oxygen saturation (St0 2 ) in the tumor regions 
[12-14]. In addition, some researchers have described tumor-to -normal contrast in tissue 
scattering parameter, ascribed to morphologic changes at a cellular level as well as changes in 
size and density of subcellular organelles in tumors [15,16]. However, it has been argued that, 
hemoglobin-DOT with no a priori information, termed standalone hemoglobin-DOT, does not 
have sufficient sensitivity for early lesion detection, and might be readily corrupted by the 
breast and its periphery structures, such as the blood vessels, glands, nipple, and chest wall 
etc., due in part to the moderate tumor-to-normal hemoglobin contrast of about 1.5-2.5 [3], as 
well as the limited quantitative accuracy and spatial resolution of the inversion [1]. To 
overcome these intrinsic shortcomings, multi-modal approaches that incorporate the 
hemoglobin-DOT with other imaging techniques have been intensively investigated, where 
the feasible-region information provided by the other modalities is normally used to 
reasonably regularize the inverse problem of the hemoglobin-DOT. Thus far, advances in 
multi-modal approaches using standard imaging modality of considerably high spatial 
resolution have been made in three fronts: concurrent optical measurements with MRI 
[17,18], X-ray mammography [19], and US [20]. 

DOT with the use of exogenous near-infrared dyes, referred to as fluorescence-DOT, was 
pioneered by Dr. Chance group for in vivo reconstruction of Indocyanine Green (ICG) uptake 
in the breast [21,22]. This modality has been attracting considerable attentions for recent years 
[23,24], in part because of its successes in molecular imaging regime of small animals [25], 
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and in part due to its contrast enhancement mechanism of significantly improving detection 
sensitivity as well as its potential to simultaneously access tissue functionality and micro- 
environmental indices such as p0 2 , pH etc. [26]. Currently, ICG, with its peak excitation and 
emission wavelengths around 780 nm and 830 nm, respectively, is the only fluorescent agent 
approved by the Food and Drug Administration (FDA) for human use. Godavarty et al. [27] 
have demonstrated use of ICG-fluorescence imaging in both canine breast cancers and 
phantoms. Corlu et al. [23] have demonstrated the first three-dimensional (3-D) ICG- 
fluorescence-DOT of in vivo human breast cancer. Other groups have also reported more 
cases of in vivo fluorescence-DOT for human breast applications [28,29]. These preliminary 
studies have observed a marked fluorescence contrast between the tumorous and normal 
tissues that originates from the agent accumulation in leaky malignant tissues [30], especially 
during the late-fluorescence phase [28]. Nevertheless, the characteristics of ICG are 
suboptimal for fluorescence-DOT of human breast: the low quantum yield often leads to 
increasingly long measuring time and auto-fluorescence that is clinically difficult to suppress - 
a challenge to instrument sensitivity and to early tumor diagnosis. Accordingly, a novel 
fluorescent agent with higher quantum yield has been under investigation by Philips since 
2002 [24]. 

While fluorescence-DOT can markedly enhance detection sensitivity of breast tumor 
imaging and is reliable in tumor localization, it relies on the uptake level of the extrinsic agent 
into the tissue that is strongly subjected to physiological and physical status of the human 
body, and lacks in quantitative merits. Therefore, from specificity point of view, hemoglobin- 
DOT that objectively discloses the pathological and physiological abnormalities of in vivo 
tissues by measuring its hemodynamic patterns is more preferred for differentiation between 
benign and malignant lesions as well as for tumor staging. To cope with both the limitations 
of hemoglobin-DOT in sensitivity and fluorescence-DOT in specificity, the combination of 
both approaches is pursued, by which the quantitative accuracy of hemoglobin-DOT is 
improved using the location a priori information provided the fluorescence-DOT, referred to 
as fluorescence-guided hemoglobin-DOT. Compared with the other multi-modality methods, 
this hybrid scheme of fluorescence and hemoglobin imaging has the potential of cost- 
effectiveness of instrumentation and an enhanced capability of detecting early lesions. 

For the above aim, we herein have developed a combined TD fluorescence- and 
hemoglobin-DOT methodology, including a 32-channel time -resolved system based on the 
time-correlated single photon counting (TCSPC) technique and a fluorescence-guided 
hemoglobin image reconstruction scheme within the framework of the generalized pulse 
spectrum technique (GPST). The performance of the methodology was evaluated using 
phantom experiments. We first assessed the linearity and sensitivity of the system responses 
to optical and fluorescence inclusions of various concentrations using two-dimensional (2-D) 
scheme, respectively. Then, the feasibility of 3-D fluorescence- and hemoglobin-DOT was 
justified using a realistic layered 2-D measuring and full 3-D reconstructing strategy. Finally, 
we performed a fluorescence-guided hemoglobin-DOT experiment and demonstrated that the 
reconstruction accuracy in hemoglobin-DOT could be significantly improved by the 
regularization of a priori fluorescence location. 

2. Methods 

2.1. Instrumentation 

Figure 1 illustrates the schematic of the developed hybrid TD fluorescence- and hemoglobin- 
DOT system. In this setup, a picosecond (ps) pulsed laser system with a multi-channel 
controller (PDL-828"Sepia II", PicoQuant, Germany) driving 2 fiber-coupled laser diode 
(LD) heads (LDH-P-780 and LDH-P-830, PicoQuant, Germany) at wavelengths of 
A, = 780 nm and A 2 = 830 nm is employed as the source. The two LD-wavelengths are 
chosen according to the peak excitation and emission wavelengths of ICG agent, respectively, 
as well as to be within the optical spectral range optimized for quantifying the tissue 
hemoglobin indexes. The controller is programmed to drive the two LD-heads in a temporally 
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Collimators Filter PMTs 
Wheels 



Fig. 1. Schematic diagram of the developed hybrid TD fluorescence- and hemoglobin-DOT 
system for breast tumor diagnosis. The system uses two /).s-pulsed LDs at 780 nm and 830 nm 
driven in a temporally multiplexing mode with the interval of 25 ns for the light delivery, and 
works with 4 PMT-TCSPC channels in a tandem serial-to-parallel scheme for the time-resolved 
detection. 

multiplexing mode with the interval of 25 ns, delivering a sequence of the alternatively pulsed 
\ - and X 1 -wavelength light with a repetition rate of 20 MHz, a power of >5 mW, and a pulse 

width of <70 ps for each wavelength. The two LD-beams are combined using a custom-made 
wave-division-multiplexer (OZ Optics, Canada) and delivered to the input fiber of a 
programmable 1 x 32-fiber-optic-switch (FSW1 x 32-MM-C, Gulin Institute of Optical 
Communications, China) that then sequentially directs the light to each of S (=32) source- 
fibers with a core diameter of 62.5 |im and a numerical aperture (NA) of 0.37. For the signal 
collection, D (=32) detection-fibers of 1000 |im in core diameter and NA = 0.22, are 
connected to four 8 x 1 -fiber-optic-switches (FSW1 x 8-MM-S, Gulin Institute of Optical 
Communications, China), whose 4 outputs are collimated and then routed to 4 PMT photon 
counting heads (PMC-100, Becker & Hickl, Germany) coupled to a 4-channel TCSPC module 
(SPC-134, Becker & Hickl, Germany) to acquire the 2-wavelength outward time-resolved 
flux, i.e. the temporal point spread function (TPSF), at 4 detector- fiber locations in parallel. 
Prior to each PMT detector is mounted a 6-hole motorized filter-wheel (FW102B, Thorlabs, 
USA) that houses 4 neutral density (ND) filters with the optical density (OD) being 1.0 to 4.0, 
respectively, and a stack of a band-pass (BP) interference filter (ICG-A Emitter, Semrock, 
USA) with transmission of >93% over the bandwidth from 813.5 nm to 850.5 nm and the 
rejection ratio of OD = 6.0 at 20 nm away from the central wavelength of 8.32 nm and a long- 
pass (LP) absorption-layered interference filter (GCC-202006, Daheng Optics, China) with 
the transmission of >85% at X = 820 nm and the rejection ratio of OD > 4.0 at X < 785 nm. 
The physical temporal-resolution of the whole system is about 250-320 ps, primarily 
depending on the transient time spread of the PMT heads that is about 180 ps and the pulse 
width of the LDs. The output-ends of the 32 source-fibers and the input-ends of the 32 
detection-fibers are paired with a side-by-side alignment to construct 32 bifurcated-fiber 
optodes. 

As depicted in Fig. 2, a cylindrical chamber with an inner diameter of 100 mm and an 
inner height of 80 mm was fabricated, within which the breast is freely suspended in the 
matching fluid. For the optode-settings, 64 holes were drilled perpendicular to the wall of the 
chamber at 4 imaging planes (named Imaging Plane #l-#4) that are evenly distributed along 
the height at an interval of 16 mm, with the 16 holes at each plane equally-spaced. The inner 
walls of the chambers are blackened to suppress the light reflection. 

2.2. Measurement protocols 

The fluorescence emission and excitation signals are detected with and without the 
fluorescence filter-stacks in place, respectively. For the TCSPC module, the time-to-amplitude 
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140mm 
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Fig. 2. Schematic diagram and photography of the cylindrical chamber (inner height: 80 mm, 
inner diameter: 100 mm). 64 holes were drilled perpendicular to the chamber wall at 4 imaging 
planes evenly distributed along the height, with the 16 equally-spaced at each plane. 

converter (TAC) range is set to 70 ns and the biased amplifier gain to 1 . Since the TAC output 
is resolved into 4096 time -bins by a 12-bit analog-to-digital converter, the above settings lead 
to a bin-width, i.e., sampling time interval, of At « 17.1 ps . 

Prior to the measurements, the realistic transmittances of the ND filters and fluorescence 
filter stacks in the 4 filter wheels are measured for the data calibration. For the former one a 
relative calibration method is adopted to preclude the influence of the dark counts on the 
measurement accuracy, where the ratio of the transmittance of the OD = n filter to that of the 
OD = n — 1 filter is measured for the 4 filter wheels all under the proper photon counts, and 
then the absolute transmittance of each filter calculated. The latter one uses the LD at \ - 

wavelength as the source to obtain the transmittances of the fluorescence filter stacks in the 4 
filter wheels, at the maximum ICG emission wavelength. 

The proposed ratio-formulated schemes of fluorescence- and hemoglobin-DOT, as will be 
described later, does not require the determination of the time-origins of the system channels, 
yet, it is still significantly beneficial to analyzing the sequencing of the measurements among 
the channels for an optimal setting of the TCSPC range as well as for reasonably choosing an 
informative and noise-robust transform-factor pair [31]. We measure approximately the time- 
origins of the system channels using the double-calibration procedure similar to that proposed 
by Hillman et al. [32]: (1) a source-calibration that collects the instrumental response 
functions (IRFs) by the first detector from the sth source (s = 1, 2, S), with their mean 
flight time denoted by (t} s 1 ; (2) a detector-calibration that acquires the IRFs by the dth 

detector (d = 1, 2, .. ., D) from the first source, with (t\ representing their mean flight time. 

With these two calibrating steps involving only 64 IRF measurements, the mean flight time of 
the IRFs for any source-detector pairs that is taken as the time-origin is calculated with a 

simple arithmetic: (t) sd + { t ) ld ~(0n • ^he calibration method and tool designed for 

the time -origin determination has been described by Gao et al. [31]. 

The workflow of the complete hybrid fluorescence- and hemoglobin-DOT process is 
shown in Fig. 3. Rigorously, it consists of five phases, referred to as ND-filter determination, 
dual-wavelength hemoglobin-DOT, fluorescence-DOT, extraction of the region of interest 
(ROI), and fluorescence-guided hemoglobin-DOT: 

Phase #1: a dual-wavelength DOT measurement is performed in the temporal 
multiplexing mode, as outlined in Sec. 2.1, with the matching fluid filled to adjust the ND 
filters for each of the used source-detector pairs so that the photon counting rate reaches the 
maximum below l/20of the LD repetition rate, i.e., 1 MHz. This upper limit for the photon- 
counting-rate can reliably prevents the "pile-up" effect of the TCSPC system during the 
successive hemoglobin-DOT measurement [33]. The phase ends with a record of the OD 
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Phase SI: ND-filter determination 



Reference-DOT data 



Reference-DOT data 



Phase #2: Dual-wavelength hemoglobin-DOT 



Hemoglobin-DOT data at A- 1 



Phase #3: Fluorescence-DOT 



Yield Image 



Images of the optical properties at A 



Phase #4: Extraction of the region of interest (ROI) 



Phase 5: Fluorescence-guided hemoglobin-DOT 



Fig. 3. Workflow of a complete hybrid fluorescence- and hemoglobin-DOT process, where the 
supporting information required by a phase is explained beside the associated arrows. 



indices of the ND filters, and a set of TPSFs simultaneously acquired at \ and Z, for all the 
source-detector pairs, referred to as the reference-DOT data set. 

Phase #2: With the subject in place and the ND filters fixed, a hemoglobin-DOT 
measurement is performed in the temporal multiplexing mode, to simultaneously acquire the 
TPSF data sets (hemoglobin data set) at the two wavelengths, as shown in Fig. 4(a). From the 
ratio of the hemoglobin (task) data set to the reference one, the background absorption and 
scattering images of the subject at the dual-wavelengths are then reconstructed for the 
successive fluorescence-DOT imaging. 

Phase #3: the fluorescence-DOT measurement is performed in the similar way to that in 
Phase #2, except that the X 1 -LD is powered off and the fluorescence filter-stacks are in place. 
The A, -LD serves as the excitation source to the ICG targets and the fluorescence TPSF data 
set around the peak emission-wavelength of ICG ( 1 2 ) is acquired, as shown in Fig. 4(b). The 
obtained fluorescence-DOT data and the hemoglobin-DOT data at A, acquired in Phase #2 

are used to calculate the Born-ratio with considering the actual integration times, the ODs of 
the used ND filters as well as the transmittances of the fluorescence BP-filter-stacks in the 
TCSPC channels, and then the yield- and lifetime-images of the fluorescence-DOT are 
reconstructed with the support of the background optical properties at X 2 . 
Phase #4: the fluorescence-yield image is segmented to extract the ROI. 




Time [ns] Time [ns] 

(a) (b) 

Fig. 4. Normalized TPSFs acquired in parallel by Detection-fiber #5-#8 for the illumination of 
Source-fiber #1 in (a) 2-D hemoglobin-DOT experiment with the reference scenario, as 
described in Sec. 3.1, and (b) 2-D fluorescence-DOT experiment with a target ICG 
concentration of 250 nM, as described in Sec. 3.2. 
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Phase #5: the image reconstruction of the dual-wavelength DOT is re -performed with the 
regularization of the ROI priors, starting from the images of the optical properties obtained 
from the dual-wavelength hemoglobin-DOT in Phase #2. The same task-to-reference ratio 
data-set as in Phase #2 is used. 

According to the optode configuration of the cylindrical imaging chamber, the 
measurement can be performed in a 2-D way, where, for each source illumination on an 
imaging plane, two groups of the 4 detection libers on the same plane are switched 
successively to the 4 parallel TCSPC channels for the photon collection. With a TCSPC 
integration time of 7} , this configuration leads to a maximum data-acquisition time of about 

32x2x7) . The measured data set can then be organized either in a 2-D way for 
reconstructing the sliced images of the fluorescent and optical properties on the 4 imaging 
planes, or in a full 3-D way for producing the volume images within the whole chamber. The 
whole setup is placed in a dark box to shield the ambient light. 

A counting threshold strategy is adopted to further optimize the measurement 
performance, where a low limit is set to the counted photon number so that the 4-channel 
parallel TCSPC process continues until the low limit is reached in the channel with the 
minimum photon-counting-rate or the integration time spent. Then the actual integration time 
used by this 4-channel parallel detection is recorded for the later data calibration. By 
deploying the 4 optodes for each parallel detection group on the geometrically symmetric 
sites, this strategy can essentially reduce the data-acquisition time. 

2.3. Data preprocessing 

Prior to the image reconstruction of the fluorescence- and hemoglobin-DOT, two procedures 
are required for the pre-processing of the raw data. Firstly, the raw data is filtered by a 5-point 
median filter following by a 5-point average filter, and the dark count of each measurement, 
which is calculated as the mean of a TPSF section before the calibrated time origin, is 
deducted the filtered curves. Then, for fluorescence-DOT, the ratio of the Laplace- 
transformed ICG-fluorescence TPSF signal, I ICG {^ d ,^ s ,0) , to the transformed excitation 

one, I At , referred to as the Born ratio [34], is calculated as the data type 

u^,e)=®(d,s) Y^ff, (i) 

where £ s and % d denote the locations of the sth source and the dth detector, respectively; /?is 
the transform-factor; @(d,s) is the gain factor of the PMT -TCSPC channel associated with 
the sth source and the dth detector, determined by the actual integration time for hemoglobin- 
DOT T x (d,s), and the used transmittance of the ND filter T ND (d,s), as well as the actual 

integration time for the fluorescence-DOT T m (d,s) and the used transmittance of the 
fluorescence BP-filter stacks T BP (d,s) 

^ > T m (d,s)-T BP (d,s) 

For hemoglobin-DOT, the ratio of the Laplace-transformed TPSF for the hemoglobin- 
DOT, (% d ,t, s ,P) , to that for the reference-DOT, i[ Ref) (% d ,t, s ,P) , referred to as the task- 
to-reference ratio, is used as the data type 
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Since the ND filters are fixed and the same integration time is used for the reference- and 
hemoglobin-DOT measurements, the above task-to-reference ratio, analogous to the Born 
ratio in fluorescence-DOT, can substantially cancel out the influence of the IRFs of the 
TCSPC channels on the reconstruction accuracy and mitigate any inconsistencies between the 
measured data and their model-predicted noiseless counterpart at high signal-to-noise ratio 
(SNR)[31]. 

2.4. Image reconstruction algorithms 

For TD hemoglobin-DOT regime with incidence of a fiber-coupled /M-pulsed laser, as 
adopted in the proposed system, the physical process of the photon migration in tissue domain 
Q is commonly described by the time -dependent diffusion equation (DE) with an isotropic 
point source in both time and space (<5-source), and its Laplace-transformed form can be 
written as [35,36] 

[-V .*-(M ; .)V (Mjc + Z?)]^ (t£„J3) = <?, (r-O, i=\, 2 (4) 

where £ s is now located at a depth of one scattering length beneath the surface; n a (r,/l, ) is 
the wavelength-dependent absorption, and /r(r,A ( .) = l/[3// s '(r,A ; )] is the diffusion 
coefficients with //'(r,A,) being the reduced scattering coefficient [37]; <t> A (r,£ s ,/?) is the 
Laplace transform of the time-dependent photon-density ct> A (r,^ s ,t), in response to the X r 

wavelength 5-source at £ s ; c is the velocity of light in tissue. 

In TD fluorescence-DOT, migration of the excitation photons follows the similar way as 
in TD hemoglobin-DOT, with the exception that the absorption coefficient is constituted of 
the intrinsic portion contributed by the tissue chromophores, and the extrinsic one originated 
from the injected ICG-agent, i.e., ju a (r,X) = H chrom (r,/l .) + H ICG (r,2 .) , while the Laplace- 
transformed fluorescence photon density, <D /CG (r,£ s ,/?) , is governed by the following 
coupled diffusion equation [38] 

[-V • *(r, A, ) V + ( Ma (r, A, ) c + fi)] 0 1CG (r, ^ ,fi) = cO^ (r, ^ , fl ^jpML ( 5 ) 

i + pt ICG \J,A\ ) 

where y ICG {r,A i ) = t] ICG (r,A l )iu ICG {r,A i ) is referred to as the yield of the ICG-agent, and 
r] ICG {y,\ ) and t icg (r,/^) are the quantum efficiency and lifetime of the agent, respectively; 

c is the velocity of light in tissue. 

The image reconstruction in hemoglobin-DOT uses a modified GPST-based nonlinear 
scheme. In this approach, a linearized formulation with regard to recovering the intermediate 

image, /f (r,/?) = cStf (r,X) + [/,<*> (r,A i )c + /]]/ M '} k) (r,A ; .)^f)(r,A,.), is firstly given 

by the iteration [35] 

/, fa&.fi.tfKrf*)-^ ( 4. .:../>)•/, (^,/U,°U' (0) ) = 

(6) 
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where j}^ and are the absorption and reduced scattering coefficients at the kth iteration, 
respectively, with //f = /uf 1 and ju'} 0 ^ = ju'}^ denoting the mean background absorption and 
reduced scattering coefficients, and Sfi^ and Sju'}^ the perturbations to ji^ and //' ( *' ; 
I^ d & s ,P,$\/jj- k) } is the model-predicted flux of the hemoglobin-DOT at the dth 

detection site % d and for the 5th source irradiation at ^ ; <S X (r,^,/?,/^*', is the 

photon density for the ^th source irradiation at £ s ; G A ^ d ,r,P,ju a k \ju'} k ^ is the flux- 
regarding Green's function for the detector at \ d . 

Then, P[ k \r,p) is reconstructed by inverting Eq. (6) with two distinct X -wavelength 
dependent transforming-factors, P + (/I,) and P_(A i ) , to separate the perturbation quantities, 
^'(r,A) and S^ k) (r,X) 

_^J^(r^) c + A (^)]^[r,AU)]-[/<! t H^) C + i94^)]^ ) ['.A(A)] 



[A (4) -A (A)]c 
^(r,lj{pf[r,^(X)]-P«[r,X,^(l)]} 



(7) 

For TD fluorescence-DOT, a GPST-based linear scheme can be directly established with 
regard to the intermediate quantity, F ICG (r,P) = y ICG (r)/[l + pr ICG (r)] , on the assumption 
that the influence of the agent absorption on the photon migration can be neglected [38] 

U^4>/?K {\^,P) = \cG h {\ d ,r,P)<$> k (r£ s ,p)F ICG (r,p)dr (8) 

where 1^ (i^,^,/?) is the boundary flux of excitation light, i.e., the \ -wavelength flux in 

the hemoglobin-DOT, at the detection site at \ d and for the source at L, s ; ct> A (r,^,/?) is the 

\ -wavelength photon density for the source at £ s , while {i sd ,r,p) the X 2 -wavelength 

flux-regarding Green's function for the detector at \ d . 

Analogously, the yield and lifetime images of the ICG agent can be separated through 
reconstructing F lCG (r,P) with a pair of the transforming-factors 

y ' CG{ ' P + F ICG (r,p + )-p_F ICG (r,p_) 

r / x F icG{ r >P + )- F ica{ r 'P-) 
T ' CG{ ' P + F ICG (r,p + )-p_F ICG (r,p_) 

For implementation, a finite-element-method (FEM) is developed under the Robin 
boundary condition for numerically solving the forward problems defined in Eqs. (4) and (5), 
from which the outward flux is proportionally associated with the boundary photon density 
according to the Fick's law [39]. With the same mesh as in the FEM forward solver, Eqs. (6) 
and (8) are discretized into linear systems, and solved using the algebraic reconstruction 
technique (ART) with the random projection access order [40], respectively, for its modest 
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memory requirement and the attained high robustness for large-scale, ill-posed linear 
inversion. 

As aforementioned, the fluorescence-DOT can provide a significantly enhanced contrast 
of 2-4 times over the hemoglobin-DOT. It is possible to segment the high-contrast yield 
image to extract the prior information on the ROI geometry, which can then be used to 
restrain the image reconstruction of the hemoglobin-DOT from the whole -region to one 
according to ROI. This "hard prior" scheme of the fluorescence-guided hemoglobin-DOT 
dramatically reduces the total number of the unknowns to be reconstructed, making the 
inverse problem better posed [41]. 

For the hemoglobin-DOT with the "hard prior" regularization strategy, the iterative 
linearization scheme corresponding to Eq. (6) can be expressed by the FEM-discretization in 
the following matrix notation 



jW =J W T 
M W = jWpW 



(10) 



,(*+!) 



. jpW + p(*) 



where P (i) and are the column vectors that numerate the optical parameter pf* (r, /?) at 
the N global nodes and NR local ROI nodes, respectively; M (t) numerates 
/, (§, , 5, , p, £\ 4 k) ) - 1, (§„ ,5, ,fi) ■ /, (§, , 5, , P, M ( : ] , M'} B) ) for all the source-detector 

combinations; J ( *' and J (i) are the global and ROI Jacobian matrix at the kth iteration, 
respectively, with its computation found in [35,38] for hemoglobin- and fluorescence-DOT, 
respectively; T is an association matrix that maps the locally indexed ROI nodes to the global 
ones 



T = 



rM), r(r p r 2 '), 
r(r 2 ,r/), T(r 2 ,r 2 '), 



r(r„r; s )' 
T(r 2 ,r' m ) 



T{r N ,r;), T(r N ,r' 2 ), ■ ■ ■ , T (r N ,r' m ) 



where T(r n r]) = 



0, r,*r! 



(11) 



where r. and r' are the positional vectors of the rth global node and the y'th ROI node, 

respectively. Owing to the high contrast nature of the fluorescence targets, ROIs can be 
simply extracted by segmenting the reconstructed ICG-yield image with a threshold 



Y ,h =yiCG+P <T K 



(12) 



1 v^l 

where y ICG =— 2>/cc( r ,) a «d ^« 



1 N 2 

—T.[yica( r i)-yicG] are the mean and 

t ,=i 



variance of the yield image, respectively; pis an attending factor that is empirically selected at 
first and then can be automatically optimized following a criterion 



(r,.)-max(j; /CG (r,.)) + £ J,cg ( r , ■) -min(j; /CG (r, .)) 



(13) 



Since the derivative of the objective function is analytically unavailable with regard to p , 

a direct line-search strategy, such as the Powell method (i.e., the second-order interpolation 
method), can used to solve the above one-dimensional optimization issue for its computational 
simplicity and reasonable convergence-rate [42]. 
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3. Experiments 



For experimental validations of the methods, a cylindrical solid phantom with a length of L = 
100 mm and a radius of R = 40 mm was made of epoxy resin (Araldite GY257, Huntsman, 
USA) mixed with Ti0 2 particles (Ti-Pure R-902+, Dupont, USA) and NIR dye (NIRD-09, 
Liaonin Huahao-Lanfan Chem.Tech. Co., Ltd., China), as shown in Fig. 5(a). The background 

optical properties of the phantom were measured to be //^ B -*/ //'^^ ) = 0.0033 / 1.09 mm -1 

and jU^P / ju'} B ^ (A 2 ) = 0.0029/0.87 mm , using a time-resolved spectroscopic (TRS) 

technique [43]. These values approximate to those of the normal breast tissue published in 
literatures [44]. For emulating breast lesions, two cylindrical holes, referred to as Hole #1 and 
Hole #2, with the same radius of R t = R 2 = 7.5 mm and different heights of Z, = 70 mm and 

Z, = 65 mm were drilled at (r x = 15 mm, f?, = 0° ) and (r, = 20 mm, 0 2 = 135° ) along the z- 

axis, respectively. For the convenience of visualization, two sets of Cartesian coordinates are 
set with X- and X'-axes intersecting the axes of Hole #1 and Hole #2, respectively. To 
achieve a 2-D scenario, Hole #1 were fully filled with the contrasting solution, and Hole #2 
inserted with a homogeneous rod that has the same optical properties as the background, 
while, for 3-D case, the two holes were partially filled and then inserted the homogeneous 
rods to form two targets with lengths of i,' * 16 mm and L 2 *16mm at different center 

heights of Z = 43 mm and Z — 53 mm, respectively. Figure 5(b) illustrates the photograph of 
the 2-D experimental setup, where 16 optodes were placed at Imaging Plane # 4, i.e., Z = 64 
mm, and Fig. 5(c) demonstrates the 3-D optode configuration, where all the 32 optodes were 
installed on the 4 imaging planes with each holding 8 located at Z = 16, 32, 48 and 64 mm, 
respectively, as indicated in Fig. 5(a). 




(a) (c) 



Fig. 5. (a) Sketch of the cylindrical solid phantom (diameter: 80 mm, length: 100 mm) and the 
coordinates used in the experiments. The phantom contains two axially-parallel cylindrical 
holes with different lengths (65 mm and 55 mm) and the same radius (7.5 mm), and can be 
configured for either a 3-D or a 2-D experimental scenario, with the contrasting solution 
partially or fully filled, respectively. Photographs of (b) 2-D and (c) 3-D optode configurations. 
In the 2-D configuration 1 6 optodes are installed on Imaging Plane #4 (Z = 64 mm), whereas in 
the 3-D configuration 32 optodes are placed on all the 4 imaging planes, with 8 optodes evenly 
distributed on each plane. 

We first performed 2-wavelength hemoglobin-DOT and ICG fluorescence -DOT 
experiments in the 2-D geometry with emphasis on their sensitivity, and then evaluated the 
efficacy of the fluorescence-guided hemoglobin-DOT scheme for improving the quantitative 
performance of the DOT image reconstruction, using the 3-D geometry. 
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3.1. Hemoglobin-DOT 

The optical-contrasting solution was prepared from mixing Intralipid-10% and India-ink with 
their concentrations controlled for the specific absorption and scattering contrasts in this 
study, and its optical properties determined by the TRS method. The experiments were 
performed for the 2-D scenario, where the absorption coefficient of the target was varied with 

IT) I 

the values of 0.012, 0.010, 0.008, 0.006 mm while the scattering coefficient was 

kept at a constant of /v'' r ' « 2.0 mm 1 . 

After performing Phase #1 and #2 in the workflow described in Sec. 2.2, we firstly 
acquired the dual-wavelength reference-DOT measurement and then the hemoglobin-DOT 
(target) measurement, from which the ratio of the target data to the reference one finally is 
calculated. Because of the possibility of the failure of the diffusion approximation for the 
source-detector pairs with short spacing as well as of the influence of the source reflection 
from the inner wall of the chamber on the closely-located detectors, the data from the 9 
detectors that are opposite to the each illuminating source were used for the image 
reconstruction, leading to a 16 x 9-dimensional data set for each of the wavelengths. With an 
integration time of 10 s and a low counting limit of 10 4 , one acquisition takes a maximum of 
640 s, ignoring the time for the fiber switching. 

The nonlinear reconstruction process employed a 2-D mesh of 2400 triangular elements 
connected at N= 1261 nodes, and was terminated after 20 iterations after which no evident 
improvement in the measurement error (i.e., the difference between the measured data and the 
calculated one) would be further achieved. The /L, -wavelength transform-factors are 

empirically set to J3 ± (A t ) = M).25^iy (/L ( )c . To suppress the artifacts, a median filter 
operating on the adjacent nodes was employed after each updating stage. Figure 6(a) 
illustrates the 2-D absorption ( ju a )- and scattering ( /u[ )-images reconstructed at 780 nm and 

830 nm, respectively, with a series of the aforementioned target contrasts, and Fig. 6(b) shows 
their profiles along the X-axis, referred to as the X-profiles. 




X |mni] X |mm] X [mm] X [mm] X Imm] X [mm] 



(a) (b) 

Fig. 6. 2-D standalone hemoglobin-DOT reconstructions of a single target with varying target 
absorption coefficient of fi M « 0.012, 0.010, 0.008, and 0.006 mm , but a constant target 

scattering coefficient of fi ' a 2.0 mm ' : (a) The // o - and //'-images reconstructed at 780- 
nm and 830-nm wavelengths, respectively; (b) the X-profiles in the jU - and fj -images at 
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both wavelengths, respectively. The circles in the images indicate the correct position and size 
of the target. 



3.2. Fluorescence-DOT 

The performance of fluorescence-DOT was experimentally assessed with the 2-D optode 
configuration. To mimic the fluorescing targets, a mixture solution of Intralipid-1% and ICG- 
dye (Cardiogreen-12633, Sigma-Aldrich, USA) was made, with the controlled ICG 
concentrations of 250 nM, 125 nM, 62.5 nM, and 30 nM, leading to a decreasing fluorescent 



yield of y] CG = 1.2x10 4 mm 

respectively, and with a Intralipid concentration for //'' r ' (A, ) = 2.0 mm 



6.0x10 " mm 



3.0xl0" 5 mm" 1 and 1.44x10^ mm" 

l(T), 



ICG Concentration = 
250nM 125nM 62.5nM 30nM 

j"l'.jSS r!A 

■1000 40n- 

■sool "S 

I -40^- 




-250 nM 
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62.5 nM 
-30nM 



1211(1 
1000- 



Fig. 7. 2-D fluorescence-DOT reconstructions of the single target with varying ICG-dye 
concentration of 250 nM, 125 nM, 62.5 nM, and 30 nM: (a) The reconstructed y - and r - 

images; (b) the X-profiles in the y - and r -images, respectively. The circles in the images 
indicate the correct position and size of the target. 

We performed firstly the reference-DOT measurement and then the fluorescence-DOT 
one, from which the normalized Born data were calculated. The same 2-D FEM meshing 
strategy as in the above 2-D hemoglobin-DOT was used for image reconstruction. To achieve 
a reasonable balance between the measurement SNR and the data-acquisition time, an 

integration time of 7^=20 8 and a low counting limit of 5xl0 3 were set for the 
fluorescence-DOT measurement, while the same settings as the previous Subsection 3.1 was 
used for the hemoglobin-DOT measurement. Figure 7 shows the 2-D yield ( y ICG )- and 

lifetime ( t icg )-images, as well as their X-profiles reconstructed with a transform-factor pair 
of p ± =±0.25//< s) (/L,)c. 

3.3. Fluorescence guided hemoglobin-DOT 

We used the 3-D configuration to demonstrate the performance of the combined diffuse 
fluorescence and optical tomography. The inversion process employed a 3-D mesh that 
contains 19200 hexahedral elements joined at N = 11349 nodes, and adopted the same 
transform-factor pairs as those in the above 2-D assessments for the fluorescence and optical 
reconstructions, respectively. To simultaneously mimic the endogenous optical contrast due to 
tumor vasculature, and the exogenous fluorescence contrast due to ICG-administration in 
breast lesions, we prepared a hybrid fluorescence- and optical-contrasting solution by mixing 
(diluted) Intralipid- 10% with India-ink and ICG with a controlled India-ink concentration of 
-0.30 mg/100 ml and a ICG concentration of -125 nM. This recipe finally results in a target 

optical contrast of Ma/Ms -0.005/1.5 mm" 1 (at both \ - and X 2 - wavelengths) and a 

target fluorescence contrast of yjJ G - 6xl0" 5 mm -1 . 



#182189 - $15.00 USD Received 21 Dec 2012; revised 21 Jan 2013; accepted 23 Jan 2013; published 25 Jan 2013 
(C) 2013 OSA 1 February 2013 / Vol. 4, No. 2 / BIOMEDICAL OPTICS EXPRESS 344 



Z=38mm Z=48mm 




Fig. 8. 3-D fluorescence-DOT reconstruction of the dual targets from the fluorescence-guided 
hemoglobin-DOT experiment: (a) The reconstructed y - and r -images of the sectional 
slice at Z = 3 8 mm and the coronal slice at Y = 0 mm, all across the center of Target # 1 , as well 
as of the sectional slice at Z = 48 mm and the coronal slice at Y = 0 mm, all across the center of 
Target #2; (b) the X-profiles in the y - and r -images of the sectional slice at Z = 38 mm 

and the X'-profiles in the y - and r -images of the sectional slice at Z = 48 mm, 

respectively. The circles and rectangles in the images indicate the correct positions and sizes of 
the targets; (c) The ROI masks of the two sectional slices at Z = 38 mm and Z = 48 mm, as well 

as the two coronal slices at Y = 0 mm and Y = 0 mm. The ROI were extracted from the y 
image of the 3-D fluorescence-DOT reconstruction with the threshold method using a attending 
factor of p = 0.86. 

Following Phase #1 to #3 of the workflow described in Fig. 3, we performed firstly the 
fluorescence-DOT reconstruction, from which the ROIs were then extracted by segmenting 

the y ICG -image by the threshold strategy with a manually optimized /?(=0.86) from the one - 
step Powell line search. Figure 8(a) illustrates the reconstructed y ]CG - and t icg -images of the 
sectional slice at Z= 38 mm and coronal slice at Y= 0 mm, all across the center of Target #1, 
as well as of the sectional slice at Z = 48 mm and coronal slice at Y = 0 mm, all across the 
center of Target #2. To facilitate the performance evaluation, the X(X')-profiles with respect 
to the above two sectional slices are plotted vs. the original ones in Fig. 8(b). Figure 8(c) 
shows the ROI masks of the two sectional and two coronal slices across the target centers. 
Finally, the hemoglobin-DOT reconstruction was performed with and without the ROI 
regularization for a comparison, with the reconstructed /u a - and /j! -images of the two 
sectional and two coronal slices across the target centers shown in Fig. 9(a) and Fig. 9(b), 
respectively, and the X(X')-profiles with respect to the two sectional slices plotted vs. the 
original ones in Fig. 9(c). 
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Fig. 9. (a) 3-D fluorescence-guided and (b) standalone hemoglobin-DOT reconstruction of the 
dual targets from the fluorescence-guided hemoglobin-DOT experiment: The /J - and // - 
images of two sectional slices at Z = 38 mm and Z = 48 mm and as well as of two coronal 
slices at Y = 0 mm and Y = 0 mm, reconstructed at 780-nm and 830-nm wavelengths, 
respectively. The circles and rectangles in the images indicate the original positions and sizes 

of the targets; (c) The X- and X'-profiles in the /J - and fj -images of the two sectional slices 

at Z = 38 mm and Z = 48 mm, reconstructed from both the 3-D standalone and fluorescence- 
guided hemoglobin-DOT at 780-nm and 830-nm wavelengths, respectively. The profiles are 
plotted vs. their respective original ones for a comparison. 



4. Discussion 

Hemoglobin-DOT reconstruction without the support of a priori information has been 
reported to be quantitatively size- and contrast-dependent [45]. According to Fig. 6(b) in 
Subsection 3.1, the fi -images reconstructed with the 2-D standalone hemoglobin-DOT 

exhibit a decreasing quantitative accuracy (defined as a ratio of the peak value in the 
reconstructed incremental target absorption or scattering to the original one) of 52%(39%), 
34%(27%), 25%(21%) and 21%(20%) at \ ( \ )-wavelength for an increasing target 

absorption contrast of 0.006, 0.008 , 0.010 and 0.012 mm -1 , respectively, demonstrating a 
saturating trend that is consistent with the previous observations in DOT investigations. In 
breast tumor diagnosis scenarios, the vascular-abnormality in cancerous regions often 

provides an absorption contrast of about 1.5-2.5, i.e., ftp - 0.004 -0.009mm -1 , depending on 
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the malignancy stage and the wavelength. This means that the standalone hemoglobin-DOT 
method can only achieve a quantitative accuracy of about 30%-50%, in terms of the above 2- 
D absorption reconstructions. For the 3-D reconstructions, as illustrated in Fig. 9, this 
quantitative accuracy is even further degraded due to smaller target-to-background volume 
ratio. It is worth noting that, regardless of the target absorption contrast, the reconstructed 2-D 
ju' s -images achieve a nearly invariant quantitative accuracy around -40% and -30% for \ - 
and X 2 -wavelengths, respectively, indicating slight crosstalk between the two parameters. 

It is well known that the ICG lifetime is about 650 ps in aqueous solution. Our previous 
investigations have observed a slight drop of the ICG lifetime to about 550 ps in Intralipid 
solution [31]. Since the background fluorescence yield is essentially null in the used phantom, 

the assumption on the background lifetime could be arbitrary from mathematical point of 
view. For the fluorescence-DOT reconstruction here, it seems by our attempts that an 

empirical value of r - 1000 ps leads to a robust and acceptable solution. As a result, both 
the 2-D and 3-D experimental results in Fig. 7 and Fig. 8 show that, despite of evident 
artifacts observed in the images due to the uncertain target-to-background contrast, the 
reconstructed lifetime roughly localizes the target but shows smaller value than the expected. 
In contrast to the lifetime reconstruction, the reconstructed yield-image correctly discloses the 
target size and location with high sensitivity (down to 30 nM for the 2-D assessments shown 
in Fig. 7), and thus can be reliably employed for the ROI extraction. 

As qualitatively demonstrated in Fig. 9(a) and Fig. 9(b), the images reconstructed in the 
hemoglobin-DOT has been significantly enhanced by means of the fluorescence-guided 
strategy. Furthermore, it is quantified from the X-profiles in Fig. 9(c) that, with the "hard- 
prior" regularization scheme, the quantitative accuracy of the reconstructed target is on 
average increased from 24.0% to 71.5% for the ^-reconstruction and from 22.6% to 53.4% 

for the ju' s . -reconstruction, showing about 3:1 and 2.4:1 improvement respectively as 

compared to the standalone hemoglobin-DOT. Table 1 summarizes the quantitative accuracy 
of the reconstructed targets. In addition to the improvement in the quantitative accuracy, the 
ability for the proposed time-domain detection scheme to effectively distinguish between the 
absorption and scattering further enables optimal selection of a tissue optical-index that is 
derived from the hemoglobin concentrations and the reduced scattering coefficient to best 
differentiate the tissue pathological states for tumor diagnosis [7]. 



Table 1. Quantitative accuracy of the reconstructed targets: fluorescence-guided 
hemoglobin-DOT vs. standalone hemoglobin-DOT 





Target #1 (//„,//) 


Target #2 (//,,//) 




Xi = 780 nm X2 = 830 nm 


Xi = 780 nm X2 = 830 nm 


Standalone 
Fluo. guided 


(23.5%, 31.0%) (19.1%, 17.2%) 
(76.6%, 86.2%) (66.7%, 47.8%) 


(29.4%, 27.6%) (23.8%, 14.6%) 
(75.3%, 52.3%) (67.5%, 27.4%) 



With an aim at demonstrating the efficacy of fluorescence-guided hemoglobin-DOT, we 
have just explored in the study the "hard-prior" regularization scheme for improving 
hemoglobin-DOT performance, which is intuitive to implement. The notable advantage of 
using a hard prior is that the total number of unknowns is dramatically reduced, making the 
inversion better-posed and thereby significantly enhancing the reconstruction quality, as 
shown in the above experiments. However, two major drawbacks of such a scheme are that its 
stability is critically dependent on the accuracy of the ROI a priori and the background optical 
structure (i.e., the absolute optical property distributions) is needed to prevent the 
reconstruction from being biased. Although the reliable ROI a priori is normally obtainable in 
phantom experiments thanks to the extremely high target-to-background yield contrast and the 
preset background fluorescence properties, this condition might not be assured in realistic 
situations due to the incomplete clearance and non-specific binding effect of the administrated 
agent [46]. Therefore, further fundamental studies are necessary for developing a "soft-prior" 
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regularization scheme for fluorescence-guided hemoglobin-DOT that is much more robust 
and unbiased in the presence of uncertainty in prior information [40]. In order to achieve 
acquisition of the background optical structure, some strategies must be worked out to cope 
with the highly ill -posed nature of the conventional pixel-based DOT, which has been used so 
far to map only the incremental optical properties over some baseline (reference-status), as 
noted in the context. It has been shown that, with the support of multi-modal structural 
images, the so-called "coarse-grained" DOT methods that are based on either region labeling 
or shape parameterization could provide promising solutions to the issue [47,48]. 

Comparing the images reconstructed in the 2-D imaging geometry (as shown in Fig. 6 and 
Fig. 7), it seems that the spatial resolution of fluorescence-DOT, even at the highest ICG 
concentration, is somewhat inferior to that of hemoglobin-DOT. This degradation in image 
quality might be primarily due to both the intrinsic low-SNR nature of the fluorescence 
signals and the measurement inaccuracy of the gain factors &(d,s) in the Born-ratio data set. 

Nevertheless, fluorescence-DOT working with effectively-enhanced tumor-to -normal contrast 
would in principle provide reliable tumor localization for quantitative improvement of 
hemoglobin-DOT, provided the system could be well calibrated and the measurement phase 
be optimized. 

5. Conclusions 

We have presented a combined time-domain hemoglobin- and fluorescence-DOT system 
based on PMT-TCSPC technique. To balance instrumental complexity and measurement 
speed, the system employs the 780 nm / 830 nm -wavelength division multiplexing, 32x1 
fiber-switching serial excitation and the 32x4 fiber-switching serial-to-parallel data 
acquisition. In practice, three primary conclusions might be obtained from developing this 
system. Firstly, fluorescence-DOT might take advantage of exogenously enhanced contrast to 
reliably profile the suspect regions. Secondly, the hemoglobin-DOT provides useful 
information on the background optical properties of the breast to improve the fluorescence- 
DOT accuracy. Thirdly, the information on the lesion localization provided by fluorescence- 
DOT is used to guide hemoglobin-DOT for enhanced characterization and quantitative 
assessment of the lesion. The experiments on a breast-mimicking phantom have validated the 
performances of the proposed methodology, and shown that the time -domain fluorescence- 
guided hemoglobin-DOT scheme presents a cost-effective and easy-to-operate way of 
markedly extending the capability of the optical mammography within the same optical 
radiation physics. Clinical evaluations of the proposed methodology are going on and the 
results will be reported in successive papers. 
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